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ABSTRACT 


We consider the possibility that galactic halos are composed of stellar remnants 
such as neutron stars and white dwarfs. On the basis of a simple model for the 
evolution of galactic halos, we follow the history of halo matter, luminosity, and 
metal and helium abundances. We assume conventional yields for helium and the 
heavier elements. By comparing with the observational constraints, which may be 
considered as fairly conservative, we find that, for an exponentially decreasing star 
formation rate (SFR) with e-folding time r, only values between 6 x 10®^r^2 x 
10 9 years are allowed together with a very limited range of masses for the initial mass 
function (IMF). Stax formation is allowed for 2MQ&mg,8MQ if r = 2 x 10 9 years, 
and for 4M©£m£6JlfQ if r = 10®years. For r = 6 x 10 8 years, the lower and upper 
mass limits merge to ~ 5 Mq. We conclude that, even though the ^possibility of 
neutron stars as halo matter may be ruled out, that of white dwarfs may still be 
a viable hypothesis, though with very stringent constraints on allowed parameters, 
that merits further consideration. 
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I. INTRODUCTION 


There is compelling theoretical and observational evidence for believing in the 
existence of baryonic dark matter, in addition to the usual non-baryonic dark 
matter, in the universe. According to the theory of big-bang nucleosynthesis (Yang 
et al. 1984; Kurki-Suonio et al. 1989), the baryon density of the universe in units of 
the critical density, ft#, * s limited to the range 0.01<£ft#;£0.1. On the other hand, 
luminous matter such as stars can only account for ft# £0.01. Possible sites for the 
baryonic dark matter include the intergalactic medium (IGM) between galaxies and 
galactic halos surrounding galaxies. 

The most direct evidence for the existence of baryons in the IGM comes from 
the Lyman a absorption lines in high red-shift quasar spectra. These systems are 
generally interpreted as being due to intervening gas clouds of almost pure H and 
He, distributed almost randomly in space and possibly confined by the pressure of 
a relatively tenuous surrounding gas (Sargent et al. 1980). The baryonic density 
contained in the Lyman a clouds is model dependent and may be as much as 
ft c h 2 ~ 0.02. Here h is the present value of the Hubble constant in units of 
100 kms~ l Mpc~^ . 

It is known that there are massive dark halos around spiral galaxies, even 
though the nature of the halo mass is unknown. Strong evidence for the existence of 
dark, massive halos comes from the rotation curves of spiral galaxies which remain 
relatively fiat out to 2-3 Holmberg radii (see, for example, Krumm and Salpeter 
1979). The inferred halo mass is several times larger than either the luminous mass 
consistent with standard population synthesis models based on a solar neighborhood 
initial mass function or the mass required to account for the dynamical mass within 
a Holmberg radius. 

The possibility that galactic halos are composed of some form of baryonic matter 
has been previously considered. For example, Hegyi and Olive (1983 and 1986) 
discussed the halos of spiral galaxies composed of snowball-like objects, gas, and 
low-mass stars or Jupiter-like objects. According to them, snowball-like objects 
should evaporate within the collapse time scale of the halos which is much shorter 
than the age of the galaxies. Halos composed of gas should heat up quickly, again, 
in a collapse time scale and emit an excessive x-ray background. If halos are made 
up of low-mass stars or Jupiter-like objects, they should produce too much infrared 
radiation to be compatible with observations, unless the slope of the IMF xj£1.7 in 
units where x — 1.35 for the Salpeter mass function. Observations tend to indicate 
that the IMF turns over (x < 0) for masses less than ~ 0.2M© (Scalo 1986). From 
the above arguments, they concluded that such forms of baryonic dark matter in the 
halo are not very plausible. Similar arguments were applied to the halo of elliptical 
galaxies by Hegyi and Olive (1989). 
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Here, we will consider the possibility that galactic halos are comprised of stellar 
remnants such as neutron stars and white dwarfs. Although Hegyi and Olive 
(1983, 1986, and 1989) also considered this possibility, they concluded that in the 
conversion of the halo mass into stellar remnants, heavy elements would be produced 
in sufficient abundances so as to over-enrich the disk. In reaching this conclusion, 
it was assumed that star formation proceeded with a constant SFR and the range 
of stellar masses was taken to be m > 2 Mq (up to ~ 100M©). It was also assumed 
that for spiral galaxies, a fraction t of the halo material is mixed into the disk. It 
was concluded that to avoid over-contamination of the disk, e < 3 x 10 -4 . For 
elliptical galaxies, e was taken to be 1 and a limit was placed on the slope of the 
IMF. In this case, over-contamination occurred if x < 10. On this basis it was 
concluded that stellar remnants are an unlikely candidate for the dark matter of 
galactic halos. 

In this paper, we will make a more detailed analysis of the evolution of the 
halo. For a constant SFR and a standard IMF, we agree qualitatively with the 
conclusions of Hegyi and Olive (1983, 1986, and 1989). Instead of setting limits 
to e, we will vary the parameters in the SFR and IMF in an attempt to find a 
consistent set of parameters such that the disk abundances of metals and helium 
are below observational bounds. We also find that a major constraint comes from 
the lower limit on the mass-to light ratio, i.e., we cannot permit large numbers of 
main-sequence stars to be still burning in the halo. In §11, we establish the basic 
equations governing the history of halo matter, metal and helium abundances in 
halos and disks, and halo luminosity. In §111, we describe several observational 
constraints which determine the plausibility of halos made up of neutron stars and 
white dwarfs. In §IV, we present our results derived from the basic equations by 
imposing these observational constraints, and in §V, we discuss various implications. 
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II. BASIC EQUATIONS 


a) Evolution of the Galactic Halo 


We begin by defining our equations and notation which follow those of Tinsley 
(1980, 1981). The creation function, or the number density of stars born over the 
mass and time intervals (m,m + dm), (t, t + dt), in the galactic halo is assumed to 
be separable, i.e., 

C(m,t) = <f>(m)ilf(t), (2.1) 

where <f> is the IMF and if> is the SFR. The IMF is normalized so that 

rmu 

I m<t>{m)dm = 1 , 

Jmi 

where n%i and mp are the lower and upper mass limits for star formation. In terms 
of <f> and "0 , the total mass densities of stars, remnants, and gas as a function of 
time are expressed as 


( 2 . 2 ) 


M s (t) = J‘[ 1 - M* ~ t')W)dt', (2.3) 

M R (t) = j‘[Mt ~ *') ~ *(< - (2.4) 

M a (t) = M t - /*( 1 - R(t - t')\ii(t')dt’ , (2.5) 

where My is the total mass density of the halo. In what follows, we take 10 12 Af© 
as the halo mass and Rjj ~ lOOAipc as the halo radius, yielding an average mass 
density Mt = PH ~ 1.6 X 10~ 3 . Here, the mass fraction of stars ever born 
which have already left the main-sequence is given by 

rmu 

/ . . m<f>(m)dm, ( 2 . 6 ) 

where mi(t) is the main-sequence turnoff mass as a function of age, and the mass 
fraction of stars returned as gas is given by 

R(t) == f mS "(m - m R )<l>(Tn)dm, (2.7) 

where m R is the mass of the remnants and rngN * s the upper mass limit for the 
supernova explosion which may differ from mjj. If mgN < m < mu, we assume 
that black holes are formed. 

We assume that the IMF of the halo stars is not very different from that of the 
disk stars and have used the IMF of the disk field stars given in Scalo (1986). From 
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his Table 7, by assuming the age of the galaxy, To, to be 15Gyr and a constant SFR 
in the disk, the IMF has been fitted to a power-law relation 


<f>(m) oc. 


( 88m 1 - 3 , 

< 0.80m- 1 - 4 , 


for 0 < m/AT© < 0.18, 
for 0.18 < tti/Mq < 1.2, 
for 1.2 < tti/Mq < 63. 


( 2 . 8 ) 


The slope of the IMF * is defined by <f> a m“( 1+ *). Notice the turnover (x < 0) for 
very low-mass stars. For the main-sequence turnoif mass, mi(t), a formula for the 
stellar lifetimes from Bahcall and Piran (1983) has been used 


f 10 — 3.61ogm + l.O(logm) 2 , for l^m/Af©^10 2 , 
\ 6.3, for m/M©^10 2 , 


(2.9) 


which agrees very closely with the values given in Scalo (1986). The mass of 
remnants (especially neutron stars and white dwarfs), mjj, has been calculated 
with the formula 



0.15m + 0.38, 
1.4, 


for 1.0 < tti/Mq < 6.8, 
for tti/Mq > 6.8, 


( 2 . 10 ) 


taken from Iben and Renzini (1983) (for their parameter rj = 1). We consider mi, 
mu, msNi <^<1 the SFR as free parameters. In particular, we use an exponentially 
decreasing SFR 

V>(0 = V’(0)exp(-t/r), (2.11) 

and allow r and V>(0) to vary. 


b) Metal and Helium Production 


If the mass fractions of metals and helium of a dying star with mass m, which 
were once processed through stellar nucleosynthesis and are returned back to the 
halo, are m 2 /m and m y /m , the total mass fraction returned to the halo in metals 
is 

tins* 

Pz(t) = I , x m z 4>(m)dm, (2.12) 

and that in helium is _ 

Py (<) = J m ^ m y <f>{m)dm. (2.13) 

Then, one can show that the metallicity in the halo, which is defined by 


z hW = 


MAt) 

MaW 


(2.14) 


where M z {t) is the mass density of metals, is expressed approximately by 


■*«=/; 


dM Stejeet (1>) 


o M G (t') dt' 


dt' 


(2.15) 
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where 

M x , ejeei (t) = J Q Pz(t - t'W)dt' (2.16) 

is the total mass density of metals ever ejected from dying stars. In deriving the 
above relation, we have assumed that z changes little between birth and death of 
the stars and that Z y. We use the approximation that the metallicity of the 
star at the time of its birth is equal to the metallicity of the halo gas at that time 
and that these elements are returned to the halo at the death of the star. We have 
also made the assumption that R(i) 3> P*(<), which is generally true. The above 
approximation is a weaker one than the commonly-used instantaneous recycling 
approximation, and we do not expect it to produce any serious errors for most of 
the cases of interest. In the same way, one can show that the helium abundance in 
the halo, defined by 

_ M,(t) 

~ Mo(*V 

where M y (t) is the mass density of helium, is expressed approximately by 




where ygj is the initial, primordial helium abundance in the halo and 


(2.18) 


eject (0 — Jq Py{t ~ t )V*(^ )dt (2-19) 

is the total mass density of helium ever ejected from dying stars. We have assumed 
that the increase in yg is small compared with ygj- Again, we have used similar 
approximations and assumptions as in the case of the metals. 


The yield of metals, m z (m), has been calculated from the estimated yields of 
Arnett (1978). Although ~ 1.7 times smaller atm = 12M© and ~ 1.5 times larger 
at m = 16AT©, the yields of Arnett agree reasonably well with those of Woosley 
and Weaver (1986) for massive stars. m z has been interpolated from Table 4 of 
Arnett (1978) for m^M©, and assumed to be negligible for m£8AT©. Similarly, 
the yield of helium, m y (m), has been interpolated from Table 4 of Arnett (1978) for 
mJfclOM© and from Figure 2 of Audouze and Tinsley (1976) for 4Af©^m^l0M©, 
and assumed to be negligible for m^4Af©. The initial primordial helium abundance 
yg t i, has been set to be 0.23. 


In order to calculate the contamination of the disk by metals and helium from 
halo stars, we assume the following: i) there is no star formation in the disk; ii) 
the halo gas moves in orbits perpendicular to the disk with typical velocity vg\ iii) 
a small fraction e of metals and helium is mixed into the disk when ejecta pass 
through the disk. Then the metallicity in the disk is calculated approximately by 



(n D /n H )M T 


3 0# (cvg 

4 fl£? P-H 



( 2 . 20 ) 
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where p z g is the mass density of metals in the halo, Sg is the area of the disk, 
^d/^H the ratio of mass in the disk to that of the halo, Rg is the radius of the 
halo, and Rg is the radius of the disk. Here, £ is a factor which accounts for the effect 
of radial orbits, in addition to tangential orbits, equal to one for tangential orbits 
and much greater than one for radial orbits. If Rg ~ Rg, the above approximation 
with only tangential orbits is expected to be good and £ ~ 1. However, if Rg < Rg, 
mixing through radial orbits becomes dominant and (R^g/Rg^C^l- In the same 
way, the helium abundance in the disk is calculated approximately by 

vdW = VD,i + (sf ) - »*<)*'. (2-21) 

where ygj is the initial, primordial helium abundance in the disk. In the calcula- 
tions, we use fig /tog ~ 0.1 and vg ~ 200 km/s and set ygj equal to ygj. On 
the other hand, the disk mass ~ 10 11 Af© yields Rg ~ 20kpc, if the column density 
of the disk is eg ~ 65M©pc -2 , corresponding to the value near the Sun (Bahcall 
1984). Hence, we assume 0.1^^(J2/j/i2fl-) 2 ^l. 

An independent effect is due to the mixing efficiency of halo gas with disk 
gas. The value of the mixing efficiency e is estimated as follows: Assume that the 
halo gas is mixed into the disk by turbulent motion. Then, the mixing length is 
approximately 

i ~ ( 2 - 22 ) 

where D is the diffusion coefficient given by 

D ~ (2.23) 

Trf is the diffusion time scale, vt is the rms turbulent velocity, and At is the charac- 
teristic scale of turbulence. Using rj ~ 3 x 10 9 years, vt ~ lOOkm/s, and Af ~ 1 kpc, 
the mixing length becomes l ~ 10 kpc and the mixing efficiency becomes c ~ 0.1. 
However, the above estimate for e should be considered as a minimum value since 
there may be other mixing mechanisms such as two body interaction of particles. 
In what follows, we shall combine the preceeding arguments with an empirical pa- 
rameters, C £ (R’D/R^h) 2 ~ 0.03, but values between 0.003 ;£C C (-Rz?/I22f) 2 ^0.1 will 
also be considered. Values of C£(Rd/Rh) 2 larger than 0.1 will be so restrictive that 
we do not consider them here. 


c) Luminosity 

The evolution of the luminosity density of the halo can be expressed as 

« = j‘ <3(1 - (2.24) 

ATf%i (t) 

Q(t) = / L(m)<f>(m)dm , 

Jmi 


where 


(2.25) 
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and L{m) is the luminosity of a star with mass m. Using L(m) in the B- 
spectral bands interpolated from the tabulated mass-luminosity relation (Allen 
1976), Lg^/MT is calculated and compared with the observed mass-to-light ratio 
of the halo in the B-spectral band. Red giants will make an additional contribution 
to the halo light which we have ignored. In our model calculations, all the above 
equations have been integrated numerically using Simpson’s rule. 
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III. CONSTRAINTS 


We now describe several observational constraints which limit the possible values 
of free parameters, m^, mp, rngg, r , and The first constraint comes from the 
metallicity and helium abundances in the disk. The Sun with an age of 4.5 x 10 9 years 
has a metallicity z ~ 0.02 and a helium abundance y ~ 0.265. According to 
the empirical age- metallicity relation of the disk obtained using F dwarf stars by 
Twarog (1980), the mean metallicity of the disk does not vary much between the 
present (1.5 x 10 10 years) and 10 10 years, but decreases by about a factor of 5 between 
10 10 years and 3 x 10 9 years. On the other hand, the extragalactic observations of 
metal and helium abundances indicate that the helium abundance could increase 
almost linearly with metallicity (Maeder 1983). For our purpose, we assume that 
the upper limits on the present metallicity and helium abundance of the disk, when 
only the contributions from the halo stars are accounted for, are 

z D ,p& 0.02, (3.1) 

and 

yD,p&VD,i + °- 05 - ( 3 - 2 ) 

Since there is also some self-enrichment of metals and helium in the disk by the disk 
stars, which is expected to be substantial, the above upper limits on zp.p and VD,p 
should be considered as fairly conservative. Also, when appropriate, we will take 
zj)( 3 x 10 9 years)£0.004. 


The next constraint comes from the gas mass density in the halo. If the halo 
predominantly consists of gas, it should be hot gas with a temperature Tg ~ 10® K 
(Spitzer 1956) and will emit x-rays. The observed, diffuse, soft x-ray background 
will therefore limit the allowed amount of halo gas density. Let C be the dumpiness 
parameter defined by 


C = 


< Pj} > 
<PB> i ’ 


(3.3) 


where pg is the baryonic mass density. From Silk (1973), the constraints based on 
the observed x-ray background in the energy range &lkeV limit 


*cn^£3, 


(3.4) 


where 6 ~ h 3 and Tg = 10 6 A can be used to set a limit on the net contribution 
of unresolved gaseous galactic halos. For tog ~ 0.1 and CQg ~ 8ttGpq/3H% ~ 
900h~ 2 MQ/MT, we get an upper limit at present, 


Mr 30h‘ 


(3.5) 


If the halo gas were to consist of smaller clouds, instead of a single massive cloud, 
we would get a more stringent limit on Mq (Tarter and Silk 1974; Hegyi and Olive 
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1986). Hence, again a fairly conservative limit would be that the present Mq should 
be less than 10% of the total mass density of the halo, or 


Mq 

Mj' 


£ 0 . 1 . 


(3.6) 


The final constraint comes from the observed mass-to-light ratio of the halo. In 
our case, the most stringent constraint is at the shortest observable wavelengths, 
since we consider stars with m^lM© which evolve into neutron stars and white 
dwarfs over a time scale less than the age of the universe. In particular, we have 
used the mass-to-light ratio in the B spectral band. From Faber and Gallagher 
(1979), on the scale of binary galaxies, the mass-to-light ratio of a spiral galaxy in 
the B spectral band is 

< 3 - 7 > 

Since most of the light comes from the disk rather than the halo, a fairly conservative 
upper limit on the present mass-to-light ratio of the halo in the B spectral band 
would be ~ 100 (M©/D©)jj, or 



(3.8) 


In Hegyi and Olive (1986), the mass-to-light ratio of NGC4565 was computed 
in the I-band. To calculate M/D, the projected mass and luminosity density were 
determined and their ratio <rjf/tr£ — M/L. The projected halo mass density was 
evaluated in terms of the 21cin rotational velocity v = 253 kma~ l (Krumm and 
Salpeter 1979). It was found that 


v 2 1 _ 3.62 x 10 9 M© 
r r kpc ^ ' 


(3.9) 


Using the surface brightness data (Hegyi and Gerber 1977), the projected luminosity 
density was fit to <ri = (a/r) + b. A 2 <r upper limit on a corresponds to a 2<r lower 
limit on M/L. In the I-band, it was found that (M/D)j > 76 (M©/D©) (Hegyi 
and Olive 1986). For comparison the 2<r lower bound on M/L in the K-band is 
(M/L) k > 38(M©/D©) (Boughn, Saulson, and Seldner 1981). We have now also 
used data on NGC4565 in the R-band (Hegyi and Gerber 1977) and performed a 
similar fit to We find a 2<r lower limit 



> 113 


Me 

I-© 


(3.10) 


Because we expect the limit on M/L in the B-band to be larger still, we feel that 
the case of (3.8) is indeed conservative. (In the case of the elliptical galaxy M87, a 
similar investigation was made in the V-band (Hegyi and Olive 1989) where it was 
found that (M/L)y > 565 (M©/D©).) 
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IV. RESULTS 


Among the free parameters, V’(O) is set by fixing the present value of Mq/Mt 
for given mi, mu, mgN, and r. However, as expected from the equations in §IIa, 
V>(0) is determined mainly by r, and, for most of the cases that we have considered, 
is approximately given by 


tK°) r 

M t 


~ 3 — 5. 


(4.1) 


In what follows, we fix the value of ip(0) by assuming Mq/Mt = 0.01 at present, 
except when we limit the free parameters from the constraint Mq/Mt& 0.1. Also 
we assume that msN — m U- (But we will return to the possibility that mg n < mu 
in §V.) Then, we are left only with three parameters, mi, mu, and r, whose values 
will be limited by the three constraints described in §111, using the equations and 
relations in §11. 


From the constraint Mq/Mt£>0.1 at present and the condition that at all times 
Mq/Mt > 0, the values of mi for given mjj and r are limited. Generally, V’(O) can 
not be too large since Mq/Mt becomes less than 0 at an intermediate time (see 
equation (2.5)). Also, V>(0) can not be too small since Mq/Mt becomes larger than 
0.1 at present even though it is always larger than 0. However, if mjr is sufficiently 
large, V’(O) can have a large enough value to suppress Mq/Mt below 0.1 at present, 
while allowing Mq/Mt to always be larger than 0, because the stars formed in 
this case die quickly due to their short lifetimes (see equation (2.9)). On the other 
hand, if mi is small, a large fraction of the halo mass is trapped in low-mass stars 
and V’(O) can not have a sufficiently large value without having Mq/Mt be smaller 
than 0. Therefore, the lower bound on mi depends sensitively on and hence 
on t. In Figure 1, the line which divides the r — mi plane into the regions with 
M g /M t &0.1 and Mq/Mt^, 0.1 is plotted for mu = 8 Af®. For the larger values 
of mu, the line shifts downwards since the halo mass locked into massive stars is 
recycled quickly. For the smaller mu, it shifts upwards. However, the dependence 
on the value of mjj is weak, since the IMF, <f>(m), falls off faster than m -1 for 
m > 1.17 Mq (see equation (2.8)) and the bulk of the stellar mass is stored into 
low-mass stars rather than massive stars. In the figure, the smaller the value of r, 
the larger the lower bound on m£. 

The constraints that zd iP ^> 0.02 and yD,p^>VD,i+ 0*05 apply mainly to the values 
of mi and mu- In our model, the metal and helium abundances increase because 
of ejection from dying massive stars (m^8Af© for metals and mj£4M© for helium). 
Since the stellar lifetimes of such massive stars are short compared to the age of 
the galaxy (see equation (2.9)), they can be approximated to die as soon as they 
are born. In that case, only the number of such massive stars is important in order 
to determine zj) p and VD,p &nd. not the formation epoch. Hence, zp p and VD,p 
depend only weakly on the value of r which determines the epoch of stax formation. 
Figure 2 shows the contours which correspond to zp p = 0.02 (long- dashed line), 
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y D ,p = 0-28 (long-dashed line), and mi = mjj (solid line) on the mi — mu plane 
for (e(Ri)/R£ f) 2 = 0.1 (Figure 2a); 0.03 (Figure 2b); 0.01 (Figure 2c); and 0.003 
(Figure 2d). It also shows the contours of zj),p = 0.01 (short-dashed line) and 
VD,p = 0.255 (short-dashed line) for comparison. Here, we have used r = 10 10 years 
and yD t i — 0.23. Even though their dependence is weak, the lines corresponding 
to zj) t p = 0.02 and VD,p — 0.28 move upwards for smaller r, since the halo gas is 
consumed more rapidly in this case and has less time to be mixed into the disk. 
From these plots, we can exclude m£^4.5M© for C c { r d/ r h ) 2 = 0.1, m^^M© 
for C £ ( r d/ r h ) 2 = 0-03, t71£^5.5Mq for C s { r dI r h ) 2 = 0*01, and m£^6.5Af© for 
C £ ( r d/ r h ) 2 = 0.003. 

From Figures 1 and 2, we can see that, for instance, r should be larger than 
6 x 10 8 years if C £ (Rd/ r H ) 2 = 0.03. However, if r is too large, the present 
value of L Hi b/ m T becomes too large to be compatible with the observed mass- 
to-light ratio. The present value of Lh^r/Mt depends on the number of stars 
that survived up to the present, which, on the other hand, depends on the epoch 
of their formation, their lifetimes, and their masses, since more massive stars are 
more luminous. Because the epoch of star formation is determined by r and only 
low-mass stars have lifetimes long enough to survive up to the present, the present 
value of Lg^/MT is basically decided by the values of r and mi and depends 
weakly on the value of my. Figure 3 shows the present Le,b/^T 118 a function of 
r for mi = 2 M® and mjj = 8M© (solid line) and for mi = 4 Af© and mu = 6 M@ 
(dashed line). The curve moves downwards for larger mi, since the number of the 
surviving stars is smaller. It moves upwards for smaller m^. If r^2 x 10 9 years, the 
present value of ^e,bI^-T is always larger than 1 /100(£©/M©) b for the allowed 
values for mi and mu and we can exclude those values of r. From Figure 3, it is 
also clear that the limit on r is not very sensitive to the constraint imposed on M/L 
when ( M/L)b 100(M©/Ir©)g. For example, there would be little difference had 
we chosen ( M/L)b > 565 (Mq/Lq)b based on the elliptical galaxy M87. 

With the constraints in §111, which we consider to be fairly conservative, only 
those values between 6 x 10 8 <^t<^2 x 10 9 years are allowed, if (c(Rd/Rh) 2 = 0 03, 
with a very limited mass range for star formation. For instance, star formation 
is allowed only for 2M©^m^8Af© if r = 2 x 10 9 years, and for 4M©^m£6Af© 
if r = 10 9 years. For r = 6 x 10 8 years, the lower and upper mass limits merge 
to ~ 5 Mq. It is important to note that these requirements imply that only a 
very small portion of the IMF mass range is allowed. Though one may argue that 
this fine tuning of the IMF could conceivably arise because of the pristine initial 
conditions of the halo, after a short time the metallicity of the halo grows and 
star formation should not be all that different from that in the disk. In this case, it 
seems quite unnatural to consider only a slice of the IMF. On the other hand, values 
between 8 x 10 8 ^t^2 x 10 9 years for Cc{Rd/Rh) 2 = 0.1, 4 x 10 8 ^r^2 x 10 9 years 
for C £ ( r d/ r h) 2 = 0.01, and 2 x 10 8 ^r^2 x 10 9 years for &{Rd! r h) 2 = 0.003 are 
allowed. 
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We have also run a model with mi = 2 Mq and mg = lOOAf© with a constant 
SFR to compare with the naive calculation by Hegyi and Olive (1983, 1986). We 
find that for zjy < 0.02, c^4 x 10 -3 which is to be compared with zg < 10 -5 and 
e£3 x 10~ 4 . In addition to being more stringent, a more severe constraint comes 
from the mass-to-light ratio. In this case ( M/L)g ~ (1/4 )(Mq/Lq)b. The halo is 
more luminous than the disk. 

Figure 4 shows an example of the time evolution of Mg /Mi, Mg/ M i, and 
Mq/Mt (solid lines), yg and yg (short-dashed lines), and Lg /Mi (long-dashed 
line) for mi = 4Afo, mg — 6Af©, r = 10 9 years, ygj — yg^ = 0.23, and 
£e(Rg/ Rg ) 2 = 0.03. Gas is initially the dominant component of the halo mass, but 
soon stars and then stellar remnants take over. At present, about 99% of the total 
mass is locked inside stellar remnants, while stars contribute a negligible amount of 
mass. The gas accounts for 1% of the halo mass and is fixed at this value in all of our 
models. The helium abundances of the halo and disk increase gradually with time, 
and at present yg p ~ 0.95 and yg tP ~ 9.28. In this case, zg = zg — 0 always since 
Mg < SMq. In Figure 5, the time evolution of zg, yg, and dyg/dzg is shown for 
mi = 3 Mq and mg = QMq (solid lines), allowed if £«(#£>/ Rgy — 0.01, and for 
mi = 3 Mq and mg = 28 Mq (short-dashed lines), allowed if (e(RD/Rg) 2 = 0.003. 
dyg/dzg initially increases rapidly and later flattens. The luminosity in Figure 4 
increases and decreases with the mass density of stars, as expected. In this case, 
the halo has a peak luminosity Lg^g/Mi ~ 40(Lq/Mq)jj and is brighter than 
Lg^g/Mi ~ 10 (Lo/M®)g for ~ O.lTo. Even though the peak luminosity and the 
duration of the bright period change for different values of the free parameters, the 
integrated luminosity over the age of the galaxy does not change very much for 
most of the cases that we have calculated, and it is approximately given by 

/ = (3 _ 3.5)%?“. (4.2) 

J Mt Mq 

This corresponds to (6 — 7) x 10 63 er<7 in the B-spectral band for a galaxy of nominal 
halo mass of 10 12 Mq, or an efficiency of conversion of rest mass to radiation of 
greater than or of order 0.3%. 
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V. DISCUSSION 


In this paper, we have considered the hypothesis that galactic halos are com- 
posed of stellar remnants such as neutron stars and white dwarfs. This possibility 
has been tested by calculating the evolution of halo matter, luminosity, and met al- 
ii city and helium abundance on the basis of a simple model, and by comparing 
them with the following three constraints: i) the present disk metallicity and he- 
lium abundance contributions from halo stars should be smaller than 0.02 and 0.05, 
respectively; ii) the present mass density of the halo gas should be less than 10% 
of the total mass density; iii) the present mass-to-light ratio of the halo in the 
B-spectral band should be larger than 100(Mq/Lq)£. 

In our evolutionary model of the halo, we have assumed that the disk of the 
galaxy forms no later than the halo. The reason for this is clear. If the disk forms 
after the period of rapid star formation in the halo, the chemical abundances of the 
disk will be far too high. For example, even in the case where metals are not formed 
in the halo as in the model shown in Figure 4, after only 1 billion years following 
the initiation of star formation in the halo, the helium abundance is ~ 40%, and 
after 2 billion years it is up to ~ 60%. Clearly the disk must form before the first 
generation erf stars die. Despite the fact that the disk forms at about the same 
time as the halo, we must require star formation in the disk to be quite different. 
That is, we must choose a different SFR and IMF. Though we do not have an 
explanation for this difference, the process of halo and disk formation is understood 
poorly enough so that we can not rigorously exclude this possibility. An alternative 
hypothesis is that the disk forms out of a fresh supply of primordial material after 
halo star formation, and the associated supply of enriched gas is mostly exhausted 
to the level of a percent or less. In this case, the time delay between halo and disk 
formation could be several billion years. 

Throughout this paper, we have taken = mjj. It is of course quite possible 
that mgjf < mp in which case all stars with masses mgjf < m < mjj end up as 
black holes without the ejection of any processed material. Typically one might 
expect that mgjf~^50MQ. However, for a first generation of stars, the absence of 
heavy elements may prevent mass loss and mg jf could be quite a bit smaller. If 
n»£ > mgjf , it is not difficult to produce a halo of black holes (this is especially 
straightforward if mi is taken to be very large). Thus as in previous investigations, 
we can not exclude black holes as a halo dark matter candidate. Although one can 
perhaps choose mgjf to be small initially, after some enrichment associated with 
one or two generations of stars, it seems unlikely that mgjf < 10 M©, and we could 
not plausibly argue for a halo exclusively containing low mass black holes. 

Based on our constraints on the metallicity and helium abundance, we conclude 
that if (c(Rd/R]j)2 = 0.03, stars should have formed only with masses between 
2Af@^m^8Af© and SFR duration between 6 x 10®years^r^2 x 10 9 years in order 
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to be compatible with the above constraints. This indicates that, even though the 
possibility of neutron stars as halo matter may be ruled out, there exists a small 
window for the possibility of white dwarfs constituting a major component of the 
galactic halo. If ( £ (Rd/Rh) 2 is smaller, it is possible that more massive stars which 
end up as neutron stars have also formed. However, considering the fact that the 
above constraints are fairly conservative, more severe limits could enable us to to 
eventually exclude the possibility of white dwarfs as well as that of neutron stars 
as halo matter. 

Our models also assume that the metallicity and increase in the helium abun- 
dance are entirely due to the halo. Because this is probably unrealistic, we have 
considered the case where Z£)^0.01 and A <^0.025, i.e., we allow the disk to pro- 
cess half of its helium and heavy elements. These results are also shown in Figure 
2. Only for the smaller value of C £ {Rd/Rh) 2 * s this difference important. For 
example, when C c (Rd/Rh) 2 ~ 0-003, my is reduced from ~ 50Af© to ~ 25 Af©. 

In Figure 5, we also see that a further constraint is available from the calculation 
of dy/dz. Unless (e( R d/Rh ) 2 1S ver 7 small and a sizable mass range is allowed in 
the IMF, we find that dy/dz would be very large, J£100. When C £ {Rd/Rh) 2 — 
0.003, then dy/dz takes more typical values ~ 1 — 3. The reason for this is 
that when ( £ (Rd/Rh) 2 ^ 0.01, the IMF contains very few stars producing metals 
and dy/dz becomes large. When a more standard IMF is used, metal producing 
massive stars exist and dy/dz falls. We note that in all of the models considered 
no deuterium is produced and thus we still require big-bang nucleosynthesis. Also, 
unless m{7^14M© there is a marked lack of oxygen in the heavy elements produced. 
In most of the cases where my < 20 Af©, we expect the chemical abundances to 
differ from those observed in extreme population II stars. 

A direct observational test of the baryonic dark matter hypothesis would be 
to search for halo white dwarfs. Models by Tamanaha et al. (1989) suggest 
that proposed surveys should be capable of detecting nearby halo white dwarfs, 
provided that the halo formed no more that 6 billion years prior to the local (solar 
neighborhood) disk. Our model, however, requires the disk to have formed soon 
after the halo and any white dwarf component should therefore be detectable with 
a modest improvement in current limits. 

It is of interest to note that the recently reported (Matsumoto et al. 1988) but 
not yet confirmed submillimeter excess in the cosmic background radiation may 
require substantial energy input of an early generation of stars. For example, Lacey 
and Field (1988) point out that the nuclear energy release associated with formation 
of the observed luminous component of galaxies falls short by an order of magnitude 
or more in providing an adequate energy source to heat up the intergalactic medium 
and distort via Comptonization the cosmic blackbody photons. Dust models for the 
submillimeter excess (Bond, Carr, and Hogan 1989; Adams et al. 1989) also require 
a very substantial stellar energy input that is more readily obtained if halo dark 
matter consists of stellar remnants. For example, with a luminous galaxy density 
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of ~ O.Olh 3 Afpc -3 , equivalent to fl ~ 0.03 /i(Mt'/ 10 1 ^M©), baryonic dark halos 
forming at redshift z* yield a diffuse radiation background, even if we only count 
the rest-frame B-band, of energy density ~ 1.3/i 3 (l + z*) _1 eV/cm 3 : for comparison, 
the submillimeter excess may amount to ~ 0.05 eV/cwr. 

In summary, we have found that only if i) the disk forms no later than the halo; 
ii) the SFR in the halo is rapidly decreasing with an e-folding time of (0.6 — 2) x 
10®years; iii) the IMF mass range in the halo contained only the segment 2 — 8Jli©; 
and iv) halo material avoided substantial mixing with the disk, can dead remnants 
such as neutron stars and white dwarfs make up the dark matter in galactic halos. In 
other words, star formation in the disk is very different from that in the halo. These 
conclusions are based primarily on the constraints from the mass-to-light ratio in 
the halo (M/L)ff^l00(MQ/LQ)g and the abundances of helium (Ay£0.05) and 
metals (zp^O.02). Our arguments presume a reasonable (0.1 — 100M©) range for 
the primordial IMF: if this assumption is relaxed, one opens the Pandora’s box of 
Jupiter-like objects (see, however, constraints discussed in Hegyi and Olive 1986) or 
very massive black holes. In this mass range, only white dwarfs survive as a possible 
candidate for baryonic dark matter, provided that a very specific form is adopted for 
the IMF and SFR. If one relaxes the constraint imposed by our assumption about 
yields, black holes would be a more likely form for baryonic dark matter. We urge 
observers to take our constraints sufficiently seriously to mount a search for any 
stellar signature of relics, such as white dwarfs, that may constitute an observable 
component of a dark baryonic halo. Ultimately, the only constraints on baryonic 
dark matter will come from observations. Of course, non-baryonic dark matter (for 
a review see Primack, Seckel, and Sadoulet 1988), for which there are numerous 
particle candidates, all hitherto undetected, remains as a viable alternative. 
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FIGURE CAPTIONS 


Fig. 1. — The line which divides the r—mi plane into the regions with Mq/Mx& 0.1 
and Mq/Mt'Z, 0.1 at the present. Here, we have used mg = 8 Mq. The values 
of r and mi in the region with Mg/Mx^O.l at the present can be excluded. 

Fig. 2. — The contours which correspond to zg v = 0.02 (long-dashed line), 
VD,p = 0.28 (long-dashed line), zg p = 0.01 (short-dashed line), yD,p — 0.255 
(short-dashed line), and mi = my (solid line) on the mi — mu plane. Here, 
we have used r = 10 10 years, yg{ = 0.23, and ^e{Rg/Rg)^ = 0.1 (2a); 0.03 
(2b); 0.01 (2c); and 0.003 (2d). The region with zd )P ^ 0.02 or yg )P ?t 0.28 can 
be excluded. 

Fig. 3. — The present Lg^g/Mx as a function of r for mi = 2 Mq and mu — 8 Mq 
(solid line), and mi = 4 Mq and mu = 6 Mq (dashed line). Lgg/Mx for 
mi = 3 Mq and mu = 28 Mq is also calculated and almost identical with that 
for mi = 4 Mq and mu = 28 Mq. The present value of Lg^g/Mx is always 

larger than 1/100(-L©/M©)£ for r^2 x 10 9 years, and we can exclude these 
values of r. 

Fig. 4. — The time evolution of Mg/Mx, Mg/Mx , and Mg/Mx (solid lines), yg 
and yg (short-dashed lines), and Lgg/Mx (long-dashed line) for mi — 4 Mq, 
mu = 6 Mq, r = 10 9 years, ygj — ygj = 0.23, and C e {Rg/Rg)^ = 0.03. In 
this case, zg = zg = 0 always since Mg < 8 Mq. 

Fig. 5. — The time evolution of zg, yg , and dyg/dzg for mi = 3Mq, mg = 9 M©, 
and C e (Rg / Rh)^ — 0*01 (solid lines), and for mi = 3 M©, mg = 28 Mq, and 
(e(Rg/ Rgj 2 = 0.003 (short-dashed lines). 
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